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Context. The line profiles of the stars with Ve sin; below a few kms"' can reveal direct signatures of local velocity fields such as 
convection in stellar atmospheres. This effect is well established in cool main sequence stars, and has been detected and studied in 
' three A stars. 

' Aims. This paper reports observations of main sequence B, A and F stars (1) to identify additional stars with sufficiently low values of 

• , Ve sin ; to search for spectral line profile signatures of local velocity fields, and (2) to explore how the signatures of the local velocity 

. fields in the atmosphere depend on stellar parameters such as effective temperature and peculiarity type. 

Oh' Methods. We have carried out a spectroscopic survey of B and A stars of low Vc sin ; at high resolution. Comparison of model spectra 

with those observed allows us to detect signatures of the local velocity fields such as asymmetric excess line wing absorption, best-fit 
Ve sin i parameter values that are found to be larger for strong lines than for weak lines, and discrepancies between observed and 
modelled line profile shapes. 

Results. Symptoms of local atmospheric velocity fields are always detected through a non-zero microturbulence parameter for main 
sequence stars having T^. below about 10000 K, but not for hotter stars. Direct line profile tracers of the atmospheric velocity field are 
found in six very sharp-lined stars in addition to the three reported earlier. Direct signatures of velocity fields are found to occur in A 
stars with and without the Am chemical peculiarities, although the amplitude of the effects seems larger in Am stars. Velocity fields 
^ ' are also directly detected in spectral line profiles of A and early F supergiants, but without significant line asymmetries. 

' Conclusions. We confirm that several atmospheric velocity field signatures, particularly excess line wing absorption which is stronger 

^SJ , in the blue line wing than in the red, are detectable in the spectral lines of main sequence A stars of sufficiently low Ve sin i. We triple 

QQ . the sample of A stars known to show these effects, which are found both in Am and normal A stars. We argue that the observed line 

' distortions are probably due to convective motions reaching the atmosphere. These data still have not been satisfactorily explained by 

' models of atmospheric convection, including numerical simulations. 

' Key words, convection - stars:atmospheres - starsxhemically peculiar ~ stars:abundances ~ stars:rotation - line:profiles 

a^ ■ 
o 

J> , 1. Introduction This effect is readily visible in a large number of lower main 

' sequence stars because many such stars have the very low pro- 

/S It has been known for many years that the spectral hne profiles of jected rotational velocity values (ve sin /< 5 - 8 km s ' ) required 

^ sharp-line cool main sequence stars similar to the Sun are asym- for the effect to be detectable at reasonable resolving power 

- - metric. When the observed line profiles of such stars are com- (~ 10^) and signal-to-noise ratio (SNR ~ 200-500). A large sur- 

pared to symmetric theoretical profiles, it is apparent that the ob- vey of bisector curvature among cool stars (spectral type mid-F 

served profiles in general have a slightly deeper long-wavelength and later) of luminosity classes II-V is reported by Gray (2005]l. 
wing compared to the short-wavelength wing. This asymmetry Such line asymmetry is generally interpreted as evidence for 

IS often described by means of a line bisector (a curve bisecting t^e presence of an asymmetric convective velocity field in the 

horizontally the observed absorption profile). Such line bisec- gt^jj^j. photosphere, with most of the atmosphere rising slowly 

tors in general have the approximate form of a left parenthesis (^n^j jhus slightly blue-shifted), while over a small fraction of the 

y\ The line bisectors of these stars are usually described as photosphere there are relatively rapid downdrafts (thus slightly 

C -shaped. red-shifted). The observed combination of the spectral lines 

from these two regions leads to integrated line profiles domi- 

nated by the slowly rising material but showing weak red-shifted 

Send ojfprint requests to: J. D. Landstreet absorption wings from the downdrafts, and line bisectors with a 

* Based in part on observations obtained at the Canada-France- joushlv "C" shaped profile 
Hawaii Telescope (CFHT) which is operated by the National Research F F 

Council of Canada, the Instimt National des Sciences de I'Univers of This qualitative picture is strongly supported by 3D hydro- 

the Centre National de la Recherche Scientifique of France, and the dynamic numerical simulations of the convection in the outer 

University of Hawaii. layers of cool stars. Such calculations exhibit the expected be- 

** Based in part on observations made at Observatoire de Haute haviour of large slowly rising gas flows that are balanced by 

Provence (CNRS), France. more rapid downflows over more limited regions, and in fact the 
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predicted line profiles are in very good agreement with observed 
profiles (e.g. Dravins & Nordlund 1990 ; Asplund et al. |2000l 
AUende Prieto, Asplund, Lopez & Lambert 2002 ; Steffen et al. 
120061 Ludwig & Steffen |200^ . 

Until recently little was known about the situation farther up 
the main sequence, mainly because above an effective temper- 
ature of Te ~ 7000 K, most main sequence stars have Ve sin / 
values on the order of 10^ km s ' , which is much larger than the 
~ 1 - 10 kms"' values of cool stars. This rotational broaden- 
ing completely overwhelms any bisector curvature. In a study of 
bisector curvature of cooler stars. Gray & Nagel ( 119891 noted 
that one Am (A8V) star, HD 3883, exhibits "reversed" bisector 
curvature; that is, that the depressed line wing is on the short- 
wavelength side of the line rather than on the long-wavelength 
side. They also found such reversed bisectors in several giants 
and supergiants lying on the hot side of a line in the HR diagram 
running from about FO on the main sequence to Gl for type lb 
supergiants. Gray & Nagel argued that their data are evidence 
for a "granulation boundary" in the HR diagram, with convec- 
tive motions on the hot side of this boundary being structurally 
different from those on the cool side. 

This picture was extended by Landstreet (1998), who 
searched for A and late B stars ("tepid" stars) of very low Ve sin 
and then carefully modelled the line profiles of such stars. He 
discovered two Am stars and one normal A star that show clear 
reversed bisectors, but found that this effect is absent from late 
B stars. Landstreet further established that reversed bisectors 
appear to correlate well with the presence of measurably non- 
zero microturbulence. Consistent with the results of testing the- 
oretical model atmospheres for convective instability using the 
Schwarzschild criterion, this microturbulence dies out at about 
AO. There is thus a second "granulation boundary" in the HR 
diagram, above which granulation (and apparent photospheric 
convection) is absent. 

Landstreet also showed that the line profiles of strong spec- 
tral lines in the late A (Am) stars are very different from those 
predicted by line synthesis using the simple model of symmet- 
ric Gaussian microturbulence, and argued that such line profiles 
carry a considerable amount of information about stellar photo- 
spheric velocity fields. 

Recently the first 3D numerical simulations of convection 
in A stars have been carried out by Steffen, Freytag & Ludwig 
(12006 ) (see also Freytag & Steffen 120041 and the earlier simu- 
lations in 2D by Freytag, Steffen & Ludwig p996). These 3D 
numerical simulations show line profile distortions with approx- 
imately the right amplitude (bisector spans of 1 - 2 km s"' com- 
pared to observed values about twice that large for ~ 8000 K), 
but the depressed line wing is the long-wavelength one, as in the 
solar case. There is currently no theoretical explanation of the 
observed depressed blue line wings in A stars, or indeed any- 
where on the hot side of the granulation boundary. (Note that 
Trampedach [2004| has also carried out computations of con- 
vection in the outer envelope of an A9 star, without computing 
line profiles, but at this temperature it is quite possible that the 
star lies on the low temperature side of the granulation bound- 
ary.) 

At present, empirical information about velocity fields in 
main sequence A and B stars is very limited. For most such stars, 
the only information available is what can be deduced from the 
microturbulence velocity parameter Useful local line profile in- 
formation is available for only a handful of A and late B stars. 
Furthermore, the best available theoretical models do not predict 
the observed line asymmetries. In these circumstances, it would 
be very useful to have more stars in which line asymmetries are 



visible directly, both to study how convective phenomena vary 
with stellar parameters such as mass, and spectral peculiarity, 
and to have a larger sample of observed profiles with which to 
compare the new generation of numerical simulations now being 
produced. 

Consequently, we have carried out a survey of A and B stars 
(mostly on or near the main sequence) with the goal of identi- 
fying more stars with sufficiently small Ve sin / values that the 
atmospheric velocity field affects observed line profiles. As stars 
of really small Vg sin / are found, we have modelled them to de- 
termine what symptoms of velocity fields are contained in the 
profiles. We have also determined a number of parameters, in- 
cluding effective temperature Tg, gravity log g, abundances of a 
few elements, accurate values of Vg sin /, and (where possible) the 
microturbulence parameter ^, for the full dataset presented here. 
The goal of this survey is to provide a significantly larger, and 
reasonably homogeneous, body of observations that may be used 
to expand our knowledge of velocity fields in the atmospheres of 
A and late B stars. 

2. Observations and data reduction 

The principal difficulty in detecting and studying the effects 
of photospheric convection (or other local velocity fields) on 
line profiles in A and B stars on or near the main sequence is 
that most have rotation velocities of the order of 10^ - 3 10^ 
kms In this situation, most velocity effects on the local line 
profiles are obliterated by the rotational broadening. The sin- 
gle effect that persists even at high Vg sin / is the requirement for 
non-zero microturbulence (characterised by the microturbulence 
parameter ^), needed to force abundance determinations using 
both weak and strong lines to yield the same abundance values. 
This parameter is widely thought to be a proxy for the presence 
of photospheric convection; Landstreet (1998) has shown that 
among main sequence A and B stars, ^ goes to zero at approxi- 
mately the same effective temperature (around 10500 K) where 
the Schwarzschild instability criterion applied to Kurucz (e.g. 
1979) model atmospheres shows that the atmospheres become 
marginally or fully stable against convection. (Note that convec- 
tion becomes energetically unimportant in the construction of 
such model atmospheres, because it is so inefficient that it car- 
ries negligible heat, at a substantially lower Tg value, ~ 8250 K, 
and is usually neglected above this temperature.) The microtur- 
bulence parameter is essentially the only kind of observational 
information available at present for the study of atmospheric 
convection in most main sequence A and B stars. 

However, a few tepid stars do have small enough values of 
VgSin/, less than a few kms"', for the effects of the velocity 
field of surface convection to be directly detectable in the ob- 
served line profiles. Three A stars (two Am stars and one normal 
one) in which the effects of local velocity fields are visible in 
the lines have been identified and studied by Landstreet ( 1998); 
three other late B or early A stars (two HgMn stars and one nor- 
mal star) were also found, in which neither microturbulence nor 
line profile distortions (relative to simple models) reveal any hint 
of photospheric convection. These data suggest that convection 
becomes detectable in line profiles of main sequence stars below 
Te ~ 10000 K, but clearly with a sample so small, conclusions 
must be rather provisional. An essential first step to drawing 
more secure conclusions from observations is to obtain a larger 
sample of stars in which velocity fields can be studied directly 
in the line profiles. 

In order to enlarge the sample of main sequence A and B 
stars of very low Vg sin / which may be studied, we have car- 
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ried out a substantial observational survey. In spite of the large 
number of Vg sin / values for A and B main sequence stars in the 
current literature (see for example Abt & Morrell 1995; Levato 
et al. |T996[ Wolff & Simon ITWTI Abt, Levato & Grosso [20021 
and references therein), a survey was necessary because most of 
the available observations of rotational velocities of early-type 
stars have been carried out with a spectral resolving power R of 
the order of 1-3 10"*. This resolving power is insufficient to de- 
termine the value of sin / if it is less than 10-20 km s"' . Since 
Landstreet ( 1 19981 1 found that direct detection of velocity fields in 
line profiles requires Vg sin / to be less than about 5 to 10 km s"' , 
new observations were required to identify stars satisfying this 
criterion. 

To have adequate resolving power to identify stars with 
VeSin/ values below a few kms"', we have used the high- 
resolution spectrograph Aurelie on the 1.52-m telescope at the 
Observatoire de Haute Provence (OHP) in France, and Gecko at 
the 3.6-m Canada-France-Hawaii telescope (CFHT). Both have 
R V 1.2 IQ^. Spectra obtained with these instruments can mea- 
sure Ve sin / values down to about 2 km s"^ This resolving power 
is sufficient to (barely) resolve the thermal line width of Fe 
peak elements, which have mean thermal velocities of about 
1.7 kms"', and thus FWHM line widths of about 3 kms"'. 
Both spectrographs are cross-dispersed echelle instruments, but 
in each, only one order is observed, so that observations record 
only a single wavelength window, typically 30 to 80 A in length, 
depending on the detector used. 

Observing missions were carried out at OHP in 1993, 1994 
and 1995. For all missions, a small window covering 4615 

- 4645 A was used, isolated by a custom interference filter 
Observations were carried out at CFHT in 1993 and twice in 
1995, using the same observing window, and also a window cov- 
ering 5150 - 5185 A. The data from these runs have been dis- 
cussed by Landstreet (|19981 l. Further missions were carried out 
in 2000 and 2001 at CFHT, both to look for more very sharp- 
line A and B stars, and to increase the available data for stars 
already identified as particularly interesting. For these runs, a 
larger CCD was available, and windows at 4525 - 4590 A (in 
2000) and 4535 - 4600 A (in 2001) were isolated by a grism. 

Candidate stars for observation were largely chosen from 
stars which are reported to have values of Ve sin / below about 15 

- 20 km s"' in previous surveys of rotational velocity (see refer- 
ences above). Although these surveys lack the resolving power 
to identify directly the stars of most interest here, they do con- 
stitute a first sorting which select the few percent of all early 
type main sequence stars that may be really sharp-lined; with- 
out these surveys to guide us, the required observing programme 
would have been more than ten times larger. 

Most of the stars observed are B, A, and early F main se- 
quence stars, although a few bright giants and supergiants were 
also observed. Examples from the Am and HgMn peculiarity 
groups were included in the survey. All the stars observed are 
brighter than about V ~ 7, since this is the effective limiting 
magnitude of most of the previous Vg sin / surveys available. The 
spectra obtained typically have signal-to-noise ratio (SNR) of 
the order of 200 in the continuum. Stars that are found to have 
very small Vg sin / values are likely to be members of close binary 
systems, which may not have been previously noticed, and so the 
stars of greatest interest to us, with Vg sin i below 8-10 km s ' , 
were mostly observed more than once to test for the presence of 
a secondary in the observable spectrum. 

Spectra were reduced using standard IRAF tasks: CCD 
frames were corrected for bias and bad pixels and extracted to 



ID; the ID spectra were divided by similarly extracted flat field 
exposures; a polynomial was fit to continuum points to normalise 
the continuum to 1 .0; the spectra were wavelength calibrated, us- 
ing typically 30 ThAr lines which are fit to an RMS accuracy of 
about 3 lO""* A; finally, the spectra were corrected to a heliocen- 
tric reference frame. 

The stars from this observing programme that have been 
found to have VgSin/ values below about 15 kms"', so that 
velocity field effects are potentially detectable in the line pro- 
files, are listed in Table [T] This table provides two names for 
each star, a spectral type (generally obtained through the Simbad 
database), an indication of binarity (either from the literature 
or from our own observations), Tg and \ogg (obtained by us- 
ing available Stromgren uvby and/or Geneva photometry to- 
gether with the FORTRAN calibration programmes described 
by Napiwotzki, Schonberner & Wenske (1993: uvby) and by 
Kunzli, North, Kurucz & Nicolet ( [19971 Geneva)). The final 
columns contain results from our own analysis that will be de- 
scribed below. 

Most stars in the survey were of course found to have values 
of Vg sin / large enough that we are not able to detect directly the 
local velocity fields in the rotationally broadened line profiles. 
For such stars, we nevertheless determined a few abundances, 
and accurate Vg sin ; and sometimes ^ values. Some results for 
these stars are reported in Table [T] but these stars are not dis- 
cussed individually any further in this paper. 

3. Treatment of double-line spectroscopic binaries 

Several of the stars in the present sample, noted as SB2 in 
Table [T] are observed to be double-lined spectroscopic binaries. 
To study the primaries of such systems, it is very useful to sub- 
tract the secondary spectrum from the observed composite spec- 
trum. This is not essential for determination of the rotation ve- 
locity of the primary (as long as we are careful to avoid lines 
significantly blended by the secondary star), but determination 
of the microturbulence parameter ^ rests on line-by-line abun- 
dance measurements, which in turn require that the line depths 
of the primary spectrum be correctly normalised. 

In the cases of greatest interest, the secondary is found from 
neutral/ion line ratios to be a mid or late A spectral type of very 
low Vg sin /. We have previously found (Landstreet 1 19981 1 that 
the line profiles of the stronger lines in such stars are not cor- 
rectly described by the standard line model with thermal broad- 
ening, damping, and height-independent Gaussian microturbu- 
lence. Thus, it is not practical to compute a synthetic spec- 
trum of the secondary which we could subtract from the pri- 
mary. Furthermore, for none of the systems studied here do we 
have enough information available (either from the literature or 
from our own data) to be able to determine accurately the sec- 
ondary/primary light ratio in our short observational windows. 

However, we have found that the secondary spectra in the 
four stars of particular interest are all very similar to the ob- 
served spectra of the very sharp-lined Am stars HD 108642 and 
HD 209625, previously studied by Landstreet (I1998I I. In all cases 
the secondary spectrum is almost unbroadened by rotation. This 
suggests that a simple method of correcting approximately the 
composite spectra for the effects of the secondary would be to 
Doppler shift an observed spectrum of one of these two Am stars 
to coincide with the secondary spectrum, subtract it, and renor- 
malise the resulting spectrum. Experiments have shown that this 
procedure in fact removes the secondary spectrum from the ob- 
served composite spectrum, leaving a residual noise in the differ- 
ence spectrum which is roughly proportional to the depth of the 
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Table 1. Sharp-line programme stars 



HD 


other 


spectral 


binarity 




log,? 


Ve sin ; 




log ;jsi/«H 


log )1Cr/'lH 


log «Fe/«H 






type 




(K) 


(cgs) 


(kms"') 


(kms"') 








2421 


HR 104 


A2Vs 


SB2 


10000 


4.0 


3.5 ± 1 


1.0 + 0.4 


-4.62 + 0.14 


-6.25 + 0.12 


-4.55 + 0.14 


27295 


53 Tau 


B9 HgMn 


SBl 


11900 


4.25 


4.9 ± 0.3 


0.4 + 0.4 


-4.60 + 0.11 


-5.85 + 0.11 


-5.32 + 0.14 


27962 


68Tau 


A2IVm 


SBl 


8975 


4.0 


11+0.5 


2.0 + 0.3 




-5.75 + 0.14 


-4.20 + 0.14 


46300 


13Mon 


AOIb 




9700 


2.1 


10 + 1 


4 + 0.5 




-6.80 + 0.14 


-4.90 + 0.14 


47105 


7 Gem 


AOIV 


SBl 


9150 


3.5 


11+0.4 


1.2 + 0.4 




-6.25 + 0.12 


-4.60 + 0.12 


48915 


aCMa 


AlVm 


SBl VB 


9950 


4.3 


16.5 ± 1 


2.1 +0.3 




-5.80 + 0.18 


-4.27+0.12 


61421 


aCMi 


F5IV-V 


SBl VB 


6500 


3.95 


7 + 1 


2.2 + 0.3 




-6.52 + 0.14 


-4.78+0.14 


72660 


HR 3383 


AlV 




9650 


4.05 


5.0 ±0.5 


2.3 + 0.3 


-4.20 + 0.11 


-6.09 + 0.10 


-4.35+0.12 


73666 


40 Cnc 


AlV 




9300 


3.8 


10.0 + 1 


2.7 + 0.3 




-6.15 + 0.14 


-4.55+0.14 


78316 


K Cnc 


B8 HgMn 


SBl 


13400 


3.9 


6.8 ± 0.5 


0.5 + 0.5 


-4.10 + 0.11 


-6.16 + 0.18 


-4.47 + 0.12 


82328 


eUMa 


F6IV 


SBl 


6360 


4.1 


9.0 ± 1 


2.2 + 0.2 




-6.62 + 0.14 


-4.95 + 0.14 


87737 


T] Leo 


AOIb 




9700 


2.0 


11 + 1 


5 + 0.5 




-6.65 + 0.14 


-4.80 + 0.14 


103578 


95 Leo 


A3V 


SB2 


8300 


4.0 


3 ± 1 


2.0 + 0.2 


-4.36 + 0.18 


-6.73 + 0.14 


-4.98 + 0.12 


107168 


8 Com 


A8m 




8150 


3.8 


13 ± 1 


3.0 + 0.5 


-4.45+0.14 


-5.75+0.14 


-4.15+0.14 


108642 


HR 4750 


A2m 


SBl 


8100 


4.1 


2±2 


4.0 + 0.5 


-4.57 + 0.19 


-6.30 + 0.14 


-4.55 + 0.14 


109307 


22 Com 


A4Vm 




8400 


4.1 


14.5 + 1.5 


3.2 + 0.3 




-6.50 + 0.14 


-4.63+0.14 


114330 


e Vir 


AlVs 


VB2 


9250 


3.4 


4± 1 


1.3 + 0.2 




-6.40 + 0.14 


-4.80 + 0.14 


120709 


3 CenA 


B5 HgMn 


VB 


17000 


3.7 


2 + 2 


0.8 + 0.8 


-4.10 + 0.18 


-7.10 + 0.31 


-4.20 + 0.18 


128167 


cr Boo 


F2V 




6750 


4.6 


9.7 ± 0.4 


1.6 + 0.3 




-6.65+0.18 


-4.90 + 0.14 


132955 


HR 5595 


B3V 


VB? 


16400 


4.4 


7 ± 1 


nd + nd 


-4.25+0.11 


-6.63 + 0.31 


-4.90 + 0.18 


142860 


7 Ser 


F6IV 




6250 


4.5 


11.5 ± 1 


1.5 + 0.3 




-6.55 + 0.18 


-4.80 + 0.14 


147084 


Sco 


A4II-III 




7850 


2.0 


5+0.5 


2.75 + 0.3 




-6.60 + 0.14 


-4.80 + 0.14 


149121 


28 Her 


B9.5 HgMn 




11000 


3.9 


9.2 ± 0.5 


0.5 + 0.5 




-6.16 + 0.11 


-4.48 + 0.11 


157486 


HR 6470 


AOV 


SB2 


9150 


3.65 


1.5 ± 1 


1.9 + 0.2 




-6.58 + 0.14 


-4.78 + 0.14 


160762 


L Her 


B3IV 


SBl 


17000 


3.9 


6.5 ± 0.5 


1 + 1 


-3.95+0.14 


-6.65 + 0.18 


-5.00 + 0.18 


170580 


HR6941 


B2V 


VB 


18000 


4.0 


5.5 ± 1 


nd + nd 


-4.00 + 0.14 




-4.90 + 0.14 


172910 


HR 7029 


B3V 


VB? 


18800 


4.4 


10 + 2 


nd + nd 


-4.02 + 0.14 


-6.72 + 0.18 


-4.90 + 0.31 


174179 


HR7081 


B3IVp 




17400 


3.9 


1.5 + 1.5 


1 + 1 


-4.45+0.18 


-6.58 + 0.18 


-4.93+0.14 


175640 


HR7143 


B9 HgMn 




12000 


4.0 


1.5 + 1 


0.5 + 0.5 




-5.65 + 0.12 


-4.98 + 0.11 


175687 


e Sgr 


AOII 




9400 


2.3 


9+ 1 


2.5 + 0.3 




-6.60 + 0.14 


-4.80 + 0.14 


178065 


HR 7245 


B9 HgMn 


SBl 


12300 


3.6 


2.0+ 1 


0.5 + 0.5 




-6.02 + 0.11 


-4.65 + 0.14 


178329 


HR 7258 


B3V 


SBl 


16600 


4.2 


7 + 1 


nd + nd 


-4.27 + 0.14 


-6.63 + 0.18 


-4.94 + 0.16 


181470 


HR 7338 


AOIII 


SB2 


10085 


3.92 


2.5 + 0.5 


0.5 + 0.5 




-6.62 + 0.14 


-4.82 + 0.12 


182835 


V Aql 


F2Iab 




6800 


2.0 


10 + 1 


5.0+ 1.0 




-6.50 + 0.14 


-4.80 + 0.14 


185395 


ecyg 


F4V 




6725 


4.45 


7 + 1 


2.0 + 0.3 




-6.55 + 0.12 


-4.80 + 0.12 


186122 


46 Aql 


B9 HgMn 




12900 


3.7 


1.0 + 0.5 


0.5 + 0.5 


-4.60 + 0.14 


-7.65 + 0.18 


-4.05+0.14 


189849 


15 Vul 


A4III m 


SBl 


7850 


3.7 


13 + 1.5 


4+ 1 




-6.60 + 0.22 


-4.60 + 0.22 


193452 


HR 7775 


AO HgMn 


SBl 


10600 


4.15 


1.0+ 1.0 


0.5 + 0.5 


-4.30 + 0.11 


-5.65 + 0.12 


-4.17 + 0.12 


209459 


21 Peg 


B9.5V 




10300 


3.55 


3.8 + 0.5 


0.4 + 0.4 


-4.55 + 0.11 


-6.42 + 0.11 


-4.70 + 0.11 


209625 


32 Aqr 


A5 m 


SBl 


7700 


3.75 


4.5 + 1.0 


4.0+ 1.0 
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subtracted lines, and which is generally no more than roughly 
1-2% of the continuum. This procedure also allows us to iden- 
tify with some certainty which lines of the primary spectrum are 
significantly blended with secondary lines, so that we can avoid 
including these lines in our detailed analysis of profile shapes. 

In general, it would be necessary to Doppler shift each pixel 
of the "template" secondary star individually, and then resam- 
ple the template star to the pixel spacing of the composite spec- 
trum. In our particular case, however, we have spectra of either 
HD 108642 or HD 209625 taken with exactly the same instru- 
mental settings used for the SB2 which is to be corrected, and it 
is found by experiment that it is adequate to simply shift the 
entire template spectrum by the integer number of pixels re- 
quired for the mean Doppler shift, and subtract the template from 
the composite spectrum pixel by pixel. This is possible because 
(a) the spectrum windows used are sufficiently short (less than 
100 A at 4500 A) and (b) because the pixel spacing of the data 
is small enough (about 0.016 A). Point (a) ensures that the error 
in the Doppler shift at the ends of the window is small, and point 



(b) allows a pixel-to-pixel correspondence to be found for which 
the fit of the template to the composite is very good. 

We have therefore corrected each of the SB2 spectra anal- 
ysed in detail by shifting a template spectrum by an integer num- 
ber of pixels, rescaling the intensity of the template to fit the 
secondary spectrum as well as possible, subtracting the template 
from the composite, and rescaling the residual primary spectrum 
to return the continuum level to 1 .0. An example of this proce- 
dure is shown in Figure [1] It will be seen that the correction of 
the primary spectrum for individual spectral features of the sec- 
ondary component is good but not perfect (especially near the 
red end of the spectrum); however, it is clear which lines in the 
primary are contaminated by lines of the secondary. 

It is clear that this procedure does not ensure an accurate de- 
termination of the light contribution of the secondary to the total 
composite spectrum, but only an estimate of this contribution. 
However, the uncertainty arising from this effect (essentially in 
the determination of the zero point of the flux in the extracted 
primary spectrum) is not very large. Typically it is found that 
the secondary template star contributes about 10% of the light 
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Fig. 1. An example of subtraction of the secondary in an SB2 system. The figure shows the original spectrum of HD 103578 in the 
4530 A window (middle spectrum, displaced upwards by 0. 1 for clarity); the spectrum of HD 108642 renormalised and pixel-shifted 
to match the secondary in the HD 103578 system as well as possible (top spectrum, shifted upwards from initial spectrum by 0.03); 
and the final spectrum of HD 103578A produced by subtraction of the scaled and shifted spectrum of HD 108642 followed by 
renormalisation to a continuum value of 1.0 (lower spectrum). 



to the composite. Since the lines of the template stars, which al- 
ways fit the profiles of the secondary star rather accurately, are 
virtually unbroadened by rotation and therefore quite deep (the 
Ba II line at 4554.03 A descends to below 0.2 of the contin- 
uum), the amount of light contributed by the secondary will not 
be overestimated by more than about 10% at the worst, or about 
1% of the total flux in the composite spectrum. Similarly, if we 
assume that the secondary in the composite SB2 spectrum has 
at least roughly solar composition, its strongest lines, which are 
highly saturated, should be about as deep as those of the tem- 
plate spectrum. Again, it appears that a 10% underestimate of 
the secondary's contribution to the total light is unlikely, which 
again leads to an overall uncertainty after subtraction of roughly 
1 % of the continuum in the residual spectrum. This is compa- 
rable to the uncertainty caused by correction for scattered light 
in the composite spectrum, and will not lead to a substantially 
increased error in either abundance determinations or in the mi- 
croturbulence parameter ^. 

The case of HD 1 14330 requires individual discussion. From 
speckle observations there is a visual secondary about 0.45 arc- 
sec from the primary, approximately 2 mag fainter than the pri- 
mary, and of similar or slightly higher Tg (e.g. Ten Brummelaar 
et al. I2OOOI 1. However, there is no trace of this component in 
the observed spectrum. Because image slicers were used in both 
spectrographs with which data were obtained for this project, 
light from this component almost certainly is part of the ob- 



served spectrum. Thus we suspect that the secondary spectrum 
is absent either because its lines are extremely broad, or because 
they are extremely naiTow and coincide with those of the pri- 
mary. We see no trace of broad lines in our (blue) spectra, so 
most probably the lines of the secondary are very narrow, and 
we must be concerned about how they influence our analysis. 
If weak secondary lines are present on the blue wing of each 
primary line, these contaminating lines could perhaps produce 
excess blue wing absorption in strong lines (or in most lines), 
and line profiles that do not closely match computed ones. It is 
not clear that such contamination could lead to strong lines that 
require a larger value of Ve sin / for best fit compared to weak 
lines. Clearly, caution must be exercised in interpreting spectra 
of HD 1 14330, and our conclusions cannot rest on this star alone. 



4. Modelling of spectra with ZEEMAN 

To uncover the information encoded in the observed spectral 
lines of a star about the atmospheric velocity field(s) present, 
we must compare the observed spectrum to a computed spec- 
trum which incorporates some model velocity fields, typically 
thermal broadening, the rotational velocity Vg sin /, and the stan- 
dard model of height-independent Gaussian microturbulence. 
The first step in this process, so that we can select an (approxi- 
mately) appropriate atmosphere model from a precomputed grid 
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of ATLAS 9 solar composition models (Kurucz [1979111 992] l. is 
to determine the effective temperature and gravity \ogg. 

These parameters were determined using available 
Stromgren uvbyH/3 and Geneva six-colour photometry, to- 
gether with appropriate calibrations. Temperatures determined 
using the two different photometry systems usually agree within 
roughly 1% (typically 100 K at ~ 10^ K), although of course 
the absolute values of are probably only known to within 
2-3%. The agreement of the inferred values of log ^ was worse, 
typically 0.2-0.3 dex, presumably because the Geneva system 
does not have a narrow-band filter equivalent to the Hj3 filter of 
the Stromgren system, and thus is less sensitive to gravity (for 
hot stars) and more sensitive to interstellar reddening (for cool 
stars). We generally adopted a value of log^ close to the value 
given by the Stromgren calibration. 

For SB2 stars, approximate values of and log g were found 
using the combined light. We tried to improve these values by 
varying the assumed values of and log g to see if the quality 
of the fit to our spectral windows was significantly improved by 
a somewhat different choice of these fundamental parameters, 
but for all the stars of interest we were unable to determine more 
accurately the fundamental parameters of the primary star from 
our short spectra. We therefore retained the estimates obtained 
using the combined light, which in general probably yield values 
of Te that are slightly too high, perhaps by 1-2%. 

The observed spectra of the stars in our sample have been 
modelled using the FORTRAN programme zeeman.f (version 
zabn4.f), described by Landstreet (1988) and Landstreet et al. 
( I1989I I, and compared to other magnetic line synthesis codes by 
Wade et al. (2001). This is nearly the same programme that was 
used in the previous study of line profiles of sharp-line A and 
B stars by Landstreet (11998) except for extensive modifications 
to include the effects of anomalous dispersion for a stellar at- 
mosphere permeated by a magnetic field. These modifications 
have essentially no impact on the analysis of the present (non- 
magnetic) stars. 

Zeeman.f makes extensive use of spectral line data obtained 
from the Vienna Atomic Line Database (Kupka et al. 120001 
[19991 Ryabchikova et al. |T997l Piskunov et al. 19951. 

A key feature of this programme is that it can be instructed 
to search for an optimum fit to a single chemical element in 
one (or several) observed spectra; at the same time, it finds the 
best fit radial velocity for the spectrum, and the best fit value 
of Ve sin /. This fitting is done assuming the standard model of 
height-independent Gaussian microturbulence, and the fitting se- 
quence can be carried out for a series of different values of mi- 
croturbulence parameter f to enable the best fitting f to be se- 
lected. The programme returns the value of the x^/v of the best 
fit (v is the number of degrees of freedom of the fit), as well as 
the parameters of the best fit. 

Zeeman.f can be used in several somewhat different modes. 
Two modes were mostly used in this study. In one, we optimise 
the fit to an element using all the significant lines of that element 
in one entire observed spectrum (30 to 75 A long), for a series 
of values of ^, selecting in the end the value of ^ which fits the 
spectrum best, and the corresponding abundance. This method 
only leads to a determination of ^ if the spectrum includes spec- 
tral lines having a range of strengths, at least some of which are 
saturated. If all the spectral lines are relatively weak, different 
values of ^ fit the spectrum roughly equally well, up to the value 
at which microturbulence broadens the line by more than its total 
observed width. 

In the second mode, we optimise the fit of a set of selected 
spectral lines of one element, one line at a time (i.e. in short win- 



dows, typically 1-2 A wide, around each of the selected lines), as 
a function of ^, and then plot "Blackwell diagrams" (Landstreet 
1998): inferred abundance as a function of ^ for each spectral 
line separately. If some of the spectral lines used are fairly sat- 
urated while others are weak, the resulting curves will usually 
intersect in a well-defined small region, which determines both 
the best fitting abundance of the element, and the best fit value 
of ^. The programme output using this method allows us to ex- 
amine conveniently the quality of the fit of an individual model 
spectral line to a particular observed line as a function of the 
various parameters such as Ve sin / and ^, and makes it easy to 
see whether or not the same fitting parameters are derived from 
all spectral lines of one (or several) element(s). Some Blackwell 
diagrams are shown by Landstreet ( I1998I I. 

The available high-quality spectra for each star have been 
analysed as described above. In general, we have optimised 
abundances of the principal detectable elements (some of He, 
O, Si, Ti, Cr, Fe, Ni, Zr, Ba) as a function of the microturbulence 
parameter ^, and chosen the microturbulence parameter value 
which yielded the best concordance of abundances derived from 
different spectral lines of one to three elements (usually Ti, Cr, 
and Fe, and sometimes Si) for which enough lines, with a large 
enough range of equivalent width, are available. The deduced 
value of ^ is then used to determine abundances of all other ele- 
ments. 

The best-fit values of Ve sin / and f , and the abundances of 
Si, Cr, and Fe, with our estimates of uncertainties, are reported 
in Table [1] In cases where the available range of line strengths 
in our spectra was not sufficient to allow us to determine the 
value of ^ we report "nd + nd" (for "not detected"). In a number 
of stars we were only able to obtain upper limits to the values 
of ^, usually from Blackwell diagrams, but occasionally from 
the very narrow line profiles; such cases are reported with equal 
values and uncertainties for ^ (e. g. 0.5 + 0.5 km s"'). 

Abundances are given as logarithmic number densities rela- 
tive to H, because this is essentially what we measure; for many 
stars in our sample the atmospheric abundance of He is un- 
known (and probably not close to solar), so we prefer to report 
log{nx/nu) rather than log(nx/ntot) values, which requires an ad 
hoc assumption about rius/nu. Uncertainties in the quoted abun- 
dance values are due in part to fitting uncertainties, which arise 
because abundances derived from individual spectral lines for a 
given value of ^ always show some scatter about the average, 
due to imprecise atomic data (especially gf values), to the fact 
that we do not have a correct model of strong lines with which 
to fit the observations, and to the fact that the atmosphere model 
used is only an approximate description of the actual stellar at- 
mosphere. These fitting uncertainties are typically of the order 
of 0.1 dex in size. There is an additional uncertainty from the 
fact that the fundamental parameters of each star are only known 
with a limited precision. We assume that the value of has an 
accuracy of approximately +3%, and that log g is known to about 
+0.25 dex. We have carried out a number of experiments to test 
the effect of these uncertainties on the derived abundances, by 
rederiving abundance values after changing Tg or log g. Because 
our abundances are generally obtained from the dominant ion, 
the abundances are not very sensitive to changes in either or 
log g. The uncertainties in fundamental parameters change the 
derived abundances by an amount which varies from element 
to element but is of the order of 0. 1 dex. The total uncertainty 
given for each abundance in Table [1] is then the sum in quadra- 
ture of the individual fitting uncertainty and a fixed uncertainty 
of 0.1 dex arising from fundamental parameter uncertainties. 
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Fig. 2. Microturbulence parameter ^ as a function of stellar Te 
for main sequence stars in the sample. Open symbols are normal 
stars, filled squares are Am stars, and filled upright triangles are 
HgMn stars. 



Note also that since we do not model macroturbulence, but 
include all broadening except for instrumental, thermal, and mi- 
croturbulent in our Vg sin / broadening, the values of Vg sin / for 
the coolest stars of our sample are larger than those given in 
studies which separate the rotational and macroturbulent com- 
ponents of line broadening. 

For the stars with the sharpest lines, the precise value of ^ ob- 
tained by our fitting procedure may not be the same as the value 
found by other methods, even with the same observational ma- 
terial and atomic data. This is because we fit line profiles using 
ax^ minimisation procedure rather than fitting equivalent width. 
As will be seen below (and was discovered by Landstreet l 19981) . 
for the sharp-line A stars the standard line profile model with mi- 
croturbulence does not fit the profiles of some very strong lines 
at all well, effectively showing that this is an inadequate descrip- 
tion of these lines. In this situation, it is not clear that there is 
any good reason to prefer one fitting method over the other - 
but it may well be that in some cases they yield slightly different 
values of ^. 

We have compared the basic parameters we have chosen 
for stars analysed, and the resulting abundance values, with a 
number of previous analyses available in the literature (e.g. Hill 
& Landstreet , 1993. Adelman et al. 1997; Hui- Bon-Hoa |2000l 
Pintado & Adelman 120031 Gebran et al. 120081 ). In most cases 
our value of Tg is within about 200 K, and our \ogg within 
about 0.2 dex, of values adopted by other authors. Our values 
of f are mostly in agreement to within 0.5 kms"' or less, and 
our resulting element abundances agree within +0.2 dex. These 
discrepancies are fairly typical for the diff'erences between inde- 
pendent abundance analyses of B and A stars, and probably de- 
rive from differences in fundamental parameters, detailed model 
atmospheres used, choices of g f values and wavelength regions, 
etc. The only significant systematic disagreement we find is that 
the element abundances we derive for the supergiants in our sam- 
ple are almost always of order -0.3 dex lower than the LTE abun- 
dances found by Venn (1995). The origin of this difference has 
not been identified, but does not affect the analysis of spectral 
line profiles which is the main subject of this paper 



5. Evidence of atmospheric velocity fields 

In our previous study of atmospheric velocity fields in middle 
main sequence stars (Landstreet ' 1 998), several kinds of obser- 
vational evidence of velocity fields were studied: strong lines 
whose profiles departed significantly from those predicted by 
the standard microturbulence model; significant (in some cases 
quite strong) excess absorption in the blue wings of strong 
lines; deduced rotation velocities Ve sin / which increase with 
line equivalent width; and the microturbulence parameter ^. 
Microturbulence has been known and studied for decades, but 
the other spectral line symptoms of local velocity fields were 
systematically studied in tepid main sequence stars for the first 
time in our previous study. 

One of the results of the spectrum synthesis carried out for 
the stars of this sample is a determination of the microturbu- 
lence parameter ^ (see Table [1) for all the stars of the sample 
except a few hot stars that have no saturated lines in the win- 
dow(s) observed. The derived values of ^ are plotted as a func- 
tion of stellar effective temperature in Figure |2] with separate 
symbols for normal. Am, and HgMn stars. This figure clearly 
shows the expected behaviour: ^ is zero within the uncertainties 
for Te > 10000 K, while below this temperature it rises rapidly 
towards lower to a limiting value which for most stars is ap- 
proximately 2 km s"' . 

The peculiar stars of this diagram have quite distinctive be- 
haviour In none of the Ap HgMn stars is a non-zero value of ^ 
detected, consistent with models of these stars that require ver- 
tically stratified atmospheres (e.g. Sigut |200T| ). In contrast, all 
the metallic-line Am stars have clearly non-zero values of ^. For 
the hotter Am stars (Tg > 8 500 K), the value of ^ is indistin- 
guishable from that of normal stars of the same Tf. value, but the 
cooler Am stars tend to have ^ values about 1 to 2 km s"' larger 
than the values found for normal stars. The relatively large ^ 
values found for the five cool Am stars of the present sample 
are consistent with the results of Landstreet ( 119981 ) for a smaller 
sample. 

The most striking symptom of atmospheric velocities dis- 
covered by Landstreet (1998) was the large departure of the 
line profiles of strong spectral lines in a few stars from the pro- 
files predicted by the standard model of a height-independent, 
isotropic Gaussian (microturbulent) velocity field. The velocity 
field signature takes the form of line profiles of strongly satu- 
rated lines that are much more triangular than line profiles pre- 
dicted by the standard model. Furthermore, these strong lines 
have substantial excess absorption in the blue line wing rela- 
tive to synthetic line profiles. Such pointed, asymmetric profiles 
are reminiscent of the line profiles found in slowly rotating cool 
stars, which are modelled by introducing macroturbulence. It is 
clear that such line distortions signal the presence of local at- 
mospheric motions more complex than those assumed by the 
standard microturbulent model, and provide valuable informa- 
tion about these motions. 

In the previous study, the phenomenon was only observed 
in very clear form in two stars (HD 108642 and HD 209625) 
with Te values close to 8000 K, both of which are Am stars. 
One further normal early A star (HD 72660) was found to show 
a very attenuated departure of observed line profiles from com- 
puted ones, primarily in the form of observed blue line wings 
that are slightly depressed below the computed line wings. For 
the other two stars in that study with below 10000 K, the 
Ve sin / values are too large for a significant departure from a 
purely rotational profile to be seen. All observed line profiles 
of three very sharp-line stars hotter than ~ 10000 K were 
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Table 2. Main sequence A-type stars in which evidence of velocity fields is directly detected in line profiles 
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Fig. 3. At most very weak convection signatures in spectral lines 
of stars with T^. values between 13000 and 9000 K observed in 
the 4530 - 4600 A window. From top to bottom, in order of 
decreasing Tei HD 186122 (46 Aql), HD 178065, HD 175640, 
HD 181470A, HD 214994 (o Peg, exceptionally observed only 
with R = 40000), and HD 73666. The top three stars are late 
B HgMn stars, while the lower three stars are normal A stars 
(one of which is an SB2). The solid lines show the observed 
line profiles of the strong 4571.97 A line of Ti ii (left panel) and 
two weak to moderate lines of Fe n at 4579.53 and 4580.06 A 
(right panel); the dashed lines show profiles computed with the 
standard line model of depth-independent, isotropic microturbu- 
lence, using the parameters of Table [T| 
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Fig. 4. Convection signatures (particularly the depressed blue 
line wing) in spectral lines of stars with in the range 10000 
to 7500 K, observed in the 4615 - 4645 A window. From top 
to bottom, in order of decreasing Te. HD 2421 A, HD 72660, 
HD 47105, HD 27962, HD 107168, and HD 209625. The top 
three stars are normal A stars, while the lower three stars are 
all Am stars. The solid lines show the observed line profiles 
of the strong 4634.07 A line of Cr ii (left panel) and the weak 
4637.50 A fine of Fe i and the blend of 4638.01 A of Fe i and 
4638.05 A of Fe ii (right panel); the dashed lines show profiles 
computed with the standard line model of depth-independent, 
isotropic microturbulence, using the parameters of Table[T] The 
bar in the lower left corner of the left panel represents a wave- 
length shift equivalent to 10 km s ' . 



found to be in very good agreement with the computed profiles; 
no indication of any atmospheric velocity field other than bulk 
stellar rotation is seen in the lines. As discussed in the intro- 
duction, the main goals of the present project are to find more 
examples of such line profile discrepancies, to discover whether 
the phenomenon is restricted essentially to Am stars, or is found 
more generally, and to determine more clearly over what range 
of Te directly observable velocity field symptoms are found. 

In Figures [3] through |5] we compare the result of detailed 
modelling of lines of the sharpest-lined main sequence stars 
in the present sample with observed line profiles. Two spectral 
windows are shown, reflecting the available observational ma- 
terial for the stars of the sample: data in the window around 
4560 A were obtained in 2000-01, while data in the window 
around 4630 A are from earlier observing runs. In each fig- 



ure, the line in the left panel is the strongest spectral line in 
the observed window that is essentially free from blends in the 
range of of the figure. The right panel shows a relatively 
weak line that is also largely blend-free. These observed lines 
are compared with model spectral lines computed assuming the 
standard model of velocity broadening due to the atomic ther- 
mal velocities, depth-independent, isotropic Gaussian microtur- 
bulence, and rigid-body rotation, along with appropriate instru- 
mental broadening. 

As mentioned above, in stars showing significant line asym- 
metry, it is also found that the Vg sin / value that best fits the 
strongest lines is almost always significantly larger (by up to 
several kms"') than that which best fits the weakest lines. 
Presumably this effect reflects the extra contribution of the local 
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Fig. 5. Convection signatures in spectral lines of stars with 
Je in the range 9250 to 6500 K, observed in the 4530 - 
4600 A window. From top to bottom: HD 1 14330, HD 157486A, 
HD 103578A (95 Leo A), HD 108642A, HD 185395 {9 Cyg), 
and HD 61421 (Procyon). The top star is a normal A star, the 
next two stars are normal (SB2) A stars, the fourth star is an Am 
star, and the bottom two are normal F stars. The solid lines show 
the observed line profiles of the strong 4563.76 A line of Ti ii 
(left panel) and the weak to moderate 4554.99 A line of Cr ii 
(right panel); the dashed lines show profiles computed with the 
standard line model of depth-independent, isotropic microturbu- 
lence, using the parameters of Table [1] The bar in the lower left 
corner of the left panel represents a wavelength shift equivalent 
to lOkms"'. 



macroscopic velocity field to the line width of saturated lines, 
and the true value of Vg sin / is (at most) as large as the minimum 
value found using the weakest, sharpest lines. In the figures, the 
projected rotational velocities of the synthetic lines with which 
these observed lines are compared are determined using only the 
weakest lines, and the same sin / value is used for all lines of 
each star. The value of ^ used is the value reported in Table [T] 
derived from the abundance modelling of weak and strong lines. 

These three figures clearly display the same phenomena 
identified by Landstreet (| 1998| ): no trace of macroscopic veloc- 
ity fields above about 10000 K, while below that temperature ^ 
rises rapidly from zero, and particularly sharp-lined stars show 
Ve sin / values that depend on line strength, line asymmetries, and 
significant discrepancies between computed and observed line 
profiles of strong lines. 

For each of the four stars in Figure [3] with effective temper- 
atures above about 10000 K (HD 186122, 178065, 175640, and 
18 147 OA, the first three of which are Ap HgMn stars), the value 
of ^ is not significantly different from zero. The same value of 
Ve sin / fits weak and strong lines, and in all cases the fit is very 
good to both weak and strong lines, which are clearly symmetric. 
These stars show no significant trace of local large-scale veloc- 
ity fields such as "macroturbulence". The same result was found 
for the normal star HD 209459 = 21 Peg (T^ = 10200 K) by 
Landstreet ( .19981 ). 

Of the seven stars in these figures with T^. between 9 000 K 
and 10 000 K, the four with Ve sin / below 6 km s"' (HD 2421 A, 
HD 72660, HD 1 14330, and HD 157486A) show small but clear 
deviations from the synthesized line profiles, specifically depres- 
sion of the observed blue line wing relative to the computed line 



for the strong (but not the weak) lines. In contrast, the two stars 
with Ve sin; ~ 10 kms"^ (HD 73666 and 47105), and the one 
star of VeSin/ = 7 kms"' observed with R - 40000 show no 
detectable effect, even though their ^ values are comparable to 
those of HD 2421 and HD 157486; for these stars the small ve- 
locities of the local velocity field are concealed by the rotational 
and/or instrumental broadening. 

There are six sharp-line stars in the sample with in the 
range 7500 to 9 000 K: HD 27962, HD 103578A, HD 107168, 
HD 108642, HD 189849, and HD 20 9625. (Two of these are the 
two Am stars studied by Landstreet 119981 ) In this temperature 
range the local velocity field amplitude is large enough to be 
clearly visible even with Ve sin / above 10 km s"' (HD 27962 and 
107168, but not in HD 189849 for which our observational ma- 
terial is of low quality), and, as previously found, in stars with 
sharper lines the impact of the local velocity field is really strik- 
ing. In HD 103578A, 108642 and 209625 the deviations from 
the model are clearly visible even in the weak lines used to de- 
termine Ve sin /. In these stars the discrepancy between the calcu- 
lated and observed line profiles is very substantial for the strong 
lines, both in the core shapes and in the occurrence of observed 
blue line wings that are shallower and deeper than the computed 
ones. In HD 27962 and 107168, both of which have Vg sin / val- 
ues above 10 km s~', the effects are visible but not striking. 

Finally, we also show two F stars in Figure |5] whose val- 
ues are below 7000 K (HD 185395 and HD 61421), and which 
show the well-known shape discrepancies of solar-type stars in 
the strong lines, with the red line wings of the strong observed 
lines shallower than and depressed below the model profile. Note 
that the line profiles of A-type stars with depressed blue wings 
as shown in Figure |5] are very different from the blueward hook 
phenomenon which has been observed for F-type stars such as 
Procyon and explained to originate from saturation of strong 
lines in upflow regions (Allende Prieto et al. 2002). The cool 
star blueward hook feature only affects lines close to their con- 
tinuum level and is much smaller in terms of velocities (less than 
0.5 kms"'). 

To provide a scale for the velocity amplitude required to ac- 
count for the line profile distortions, in each figure where these 
effects are apparent we include a small horizontal error bar with 
a total width of 10 km s 

A summary of all the main sequence A stars in which we 
find direct evidence of an atmospheric velocity field (beyond that 
provided by a non-zero value of ^) has been found is given in 
Table |2] Columns 3 to 6 recall the spectral type, binarity, effec- 
tive temperature and microturbulence parameter of these stars. 
In the column "depressed blue wing" the comment "weak" or 
"strong" indicates the extent to which the short-wavelength line 
wing is deeper than the long- wavelength wing. The column "line 
shape discrepancy" describes qualitatively the extent of discrep- 
ancies, other than the depressed blue line wing, between ob- 
served and calculated line profiles of the strongest spectral lines. 
The last two columns list extreme values of Vg sin / as derived 
from best fits to the weakest and strongest spectral lines in our 
spectra. As discussed qualitatively above, we see that the dif- 
ference between the best fit Ve sin / values for the weakest and 
strongest lines in our spectra seems to increase quite rapidly for 
Te below about 9000 K to a value that can be several km s"' at 
around 8000 K. (Note that HD 3883, noted by Gray & Nagel 
II1989I as having a reversed bisector, is omitted from Table |2] 
because we have no spectra of this star.) 

We notice that the direct evidence of velocity fields in the 
line profiles agrees very well with the indirect evidence provided 
by the microturbulence parameter Main sequence stars having 
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Te above about 10000 K have ^ values that are indistinguish- 
able from 0, and thus lack any evidence of small-scale velocity 
fields. These same stars have spectral lines that are in excellent 
agreement with calculated line profiles in which the broadening 
is assumed to be due only to rotation, thermal broadening, colli- 
sional damping, and instrumental resolution; these stars thus do 
not show any direct evidence of local atmospheric hydrodynamic 
motions. 

In contrast, stars with below 10000 K all have non-zero 
values of ^, and all the stars with "sufficiently" small values of 
Ve sin / also show direct line profile evidence of velocity fields. 
Local velocity fields are detectable through differences between 
profiles of observed strong spectral lines and computed lines, 
through asymmetries in the observed lines, and in most cases 
through best-fit Ve sin / values that are systematically larger for 
strong than for weak lines. Furthermore, the strongest directly 
detected symptoms of local velocity fields, observed in stars with 
Tg around 8000 K, generally coincide with the largest deduced 
values of ^. 

These results confirm the findings from a smaller dataset by 
Landstreet ( 1998 ), where it was shown that the range of effective 
temperatures for which non-zero-values of ^ are found coincides 
rather closely with the occurrence of a convectively unstable re- 
gion in the stellar atmosphere as predicted by the Schwarzschild 
criterion applied to the model atmosphere used for spectrum 
synthesis. Thus, the occurrence of convective motions in the at- 
mospheres is expected in essentially the range of T^. in which 
it is observed. The importance of the observations presented in 
this paper comes from the direct signatures of this convection in 
spectral lines, which should be able to provide significant con- 
straints on the nature of the actual flows occurring in the atmo- 
sphere. 

It may be noted that the line profile distortions relative to the 
canonical model are observed essentially only in lines of once- 
ionised atoms such as Fe, Cr and Ti, and not in the few visible 
lines of neutral atoms. This is presumably due to the fact that all 
neutral lines are relatively weak, while the effects we are inter- 
ested show themselves in strong, heavily saturated lines. 

It is of interest to compare the symptoms of local velocity 
fields found in the spectra of main sequence stars with those 
found in lines of supergiants of similar spectral types, and for 
this reason several of the stars in our survey are such luminous 
objects. It is well known that A and F supergiants require rela- 
tively large values of the microturbulence parameter (typically 
several kms"^) for successful spectrum synthesis (e. g. Venn 
[19951 Gray et al. iSOOl; Przybilla et al. |2506] ). and so we ex- 
pect to find strong evidence of velocity fields in the line profiles 
of the slowest rotators. 

Five supergiants, ranging in Tg from 9700 to 6800 K, are 
listed in Table[T]and are shown in Figure|6] As with the main se- 
quence stars, the Vg sin / values have been determined using only 
the weakest lines, and for each star the same value of Vg sin / has 
been used for all synthesized profiles. A systematic result of our 
use of only the weakest lines is that the Vg sin / values determined 
here are always significantly smaller than those found in previ- 
ous studies, where Vg sin / is fit mainly to the stronger lines. In 
the right panel of the figure an example of relatively weak lines 
shows that these profiles are fairly well fit by simple rotational 
broadening (together with thermal and instrumental broadening 
and classical microturbulence). In contrast, the observed strong 
lines of the left panel are all significantly different in shape from 
the computed lines. If we had (unphysically) allowed Vg sin / to 
be a free parameter for each spectral line separately, the strong 
lines could have been substantially better fit. In each case the 
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Fig. 6. Striking convection signatures in spectral lines of A 
and early F-type supergiants with Tg values between 9700 and 
6800 K observed in the 4530 - 4600 A window. From top to 
bottom: HD 46300 (13 Mon), HD 87737 (77 Leo), HD 175687 
(^1 Sgr), HD 147084 (o Sco), and HD 182835 (v Aql). The sohd 
lines show the observed line profiles of the strong 4588.20 A 
line of Cr 11 (left panel) and a group of weak to moderate lines, 
4595.36 and 4596.06 A of Fe i and 4595.68 and 4596.02 A of 
Fe II (right panel); the dashed lines show theoretical profiles 
computed with the standard line model of depth-independent, 
isotropic microturbulence computed with the parameters of 
Table [T] The bar in the lower left corner of the left panel rep- 
resents a wavelength shift equivalent to 10 km s 



best fit Vg sin ; value determined from strong lines is a few km s 
larger than that found from weak lines, and is similar to values 
found in other high-resolution studies. In the coolest star shown, 
HD 182835, the best fit Vg sin / value found from weak lines is 
no more than 10 km s"', while the best fit for the strongest lines 
is about 18 km s"' ! 

Thus, we find in these highly luminous stars evidence for the 
local velocity fields phenomenologically similar to that found 
in main sequence stars of similar Tg, but with larger amplitude. 
However, there is one significant difference: the supergiants do 
not show the strong line asymmetry seen in the main sequence 
stars with similar temperatures. The depression of the observed 
wings of strong lines below the computed wings is essentially 
the same on the blue and red sides of the line. This result is dif- 
ferent from that of Gray & Toner (1986), who find significant 
excess absorption in the blue wings of spectral lines of F super- 
giants. However, their data suggest that the asymmetry dimin- 
ishes towards earlier spectral types. As their earliest supergiant 
has spectral type F5 while our latest type is F2, perhaps we are 
simply seeing the asymmetry decrease below the threshold for 
easy visibility. 

A different possibility for the velocity field observed in su- 
pergiants is that the line profiles reveal the presence of pulsa- 
tions. Aerts et al. (2008) have shown that the qualitatively similar 
line profiles of an early B-type supergiant, typically modelled by 
macroturbulence, may be produced by numerous low-amplitude 
gravity-mode oscillations, rather than by large-scale convective 
motions. 

This raises the question of whether the velocity fields that we 
observe in sharp-line A-type main sequence stars might also be 
due to pulsations, but there are several reasons to doubt this. The 
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first is that the computations of Aerts et al. ( I2008I I do not show 
any significant line asymmetry, the most persistent feature we 
observed in the line profiles of A stars. Furthermore, we observe 
depressed blue line wings and other symptoms of a velocity field 
in stars significantly hotter than the hot limit of the instability re- 
gion (about 8700 K or a bit cooler; e.g. Joshi et al. l2003l l. A third 
reason for rejecting pulsations is that the asymmetry of line pro- 
files in tepid stars has a constant behaviour, reproducible through 
repeated observations of all the stars in Table |2] the excess ab- 
sorption is always in the short-wavelength line wing. It is dif- 
ficult to see why pulsation would not lead to excess absorption 
that would shift from one line wing to the other Finally, one 
of the stars in which the symptoms of a non-thermal velocity 
field are strongest, HD 209625, has been observed very care- 
fully to search for evidence of pulsations (Carrier et al. I2007I I. 
but no hint of pulsational variations was found. This is consis- 
tent with many observational studies which have concluded that 
only marginal and evolved Am stars are found to pulsate (cf. 
the review of Balona |2004| l. thus supporting theoretical results 
from diffusion theory (Turcotte et al. '2000). A few exceptional 
classical Am stars, pulsating at very low amplitude, have never- 
theless been identified (Kurtz |20001 l. In contrast, not all the nor- 
mal A stars in the centre of the instability strip are found to pul- 
sate above observational thresholds and 6 Set variables with low 
Ve sin / are found to have a very large spread in observed pulsa- 
tion amplitudes (Breger l2000l l. Thus, although we cannot com- 
pletely exclude pulsation from being present in the stars listed 
in Table |2] with Tg < 8700 K, their maximum variability level 
due to pulsation has to be very small, making pulsation a very 
unlikely candidate for large velocity fields that cause asymmet- 
rically broadened line profiles. 

It is of course possible to consider other explanations of 
the depressed line wings than atmospheric motions. The excess 
broadening is not caused by magnetic (Zeeman) broadening, as 
this would selectively broaden lines of large Zeeman sensitivity; 
instead one finds simply that all the strongest lines are broadened 
most. In particular, if the broadening were magnetic, the strong 
4634.07 A line of Cr ii would be considerably narrower than the 
nearby 4629.34 A line of Fe ii; this is not observed. Similarly, 
we can exclude the possibility that the effect is caused by hyper- 
fine splitting, because the principle isotopes of Fe, Cr, and Ti, 
in which strong line distortion is observed, are overwhelmingly 
even-even isotopes with no net nuclear magnetic moment, and 
also because the observed effect depends strongly on stellar Tg. 
We can exclude the possibility that the selective broadening and 
distortion is caused by rapid rotation seen pole on, as this would 
broaden both strong and weak lines by a similar amount. We 
note that also a direct interaction between the two members of 
each of the observed binaries in the sample is unlikely to explain 
the asymmetric profiles due to their constant behaviour over dif- 
ferent rotational periods and the varied nature of their binarity. 
Finally, it seems unlikely that circumstellar material (which can 
have quite large velocities) plays any role, as there is no observed 
emission at Ho- and H/3, such as one finds in Herbig AeBe stars. 

We conclude that the atmospheric convective velocity field 
is still the most probable explanation of the line profile peculiar- 
ities reported here. 

6. Conclusions 

The results described in the preceding sections lead to a number 
of conclusions. 

Local atmospheric velocity fields are confidently detected 
in main sequence stars having below 10000 K, but not in 



stars with higher effective temperatures. These velocity fields 
are always detected via clearly non-zero values of the micro- 
turbulence parameter ^, which is indistinguishable from zero for 
Tg > 1 0000 K, but rises quickly to a value near 2 km s" ' for lower 
Tg. A larger value of ^ is found essentially only for Am stars with 
Te between 7500 and 9000 K. 

In all stars with sufficiently sharp lines (vg sin / below about 
6 km s^'for stars with 9 000 < < 10 000 K, and below about 
12 or 13 km s"' for stars with 7 500 < < 9 000 K), there 
is a perfect correspondence between non-zero ^ and other ex- 
pressions of a local velocity field. Sharp-line stars with non-zero 
micro turbulence also exhibit some or all of (1) best-fit values 
of Ve sin / that increase with the equivalent width of the line, (2) 
significant line asymmetry, with a deeper depression of the blue 
line wing for stars with larger than about 7000 K, and a deeper 
red wing for cooler stars, and (3) profiles of strong lines that are 
clearly not correctly modelled by the usual LTE line formation 
theory if one assumes that the only line broadening mechanisms 
are instrumental broadening, thermal and pressure broadening, 
height-independent isotropic Gaussian microturbulence, and ro- 
tation with v'e sin / derived from the weakest available lines. 

Our search for more stars which clearly exhibit characteris- 
tics (2) and (3) has yielded only a few new examples (HD 2421 A, 
HD 27962, HD 103578A, HD 107168, HD 157486A, and prob- 
ably HD 1 14330). However, even this small number has tripled 
the available sample. The line profile characteristics of these 
stars are consistent with those of the very sharp-line stars already 
studied (HD 72660, HD 108642 A, and HD 209625; Landstreet 
1998), both in form and in variation with effective temperature. 
From the newly detected stars we learn that the modifications 
to spectral lines produced by local velocity fields are found in 
both Am and non-Am stars of similar Tg, although (consistent 
with the larger ^ values found in Am stars) the amplitude of line 
profile distortions seems larger in Am than in normal A stars. 
The effects directly observable in really sharp-line A-star spec- 
tra become stronger as decreases below 10000 K, down to 
about 7 000 K, below which the line profiles take on the already 
well-known shape characteristic of solar-type stars. 

There are now sufficient data available to make useful com- 
parison with numerical simulations of convection for an inter- 
esting range of values of and log g, and the data for different 
tepid main sequence stars show a reassuring degree of consis- 
tency. We should emphasize that even if the effects we are study- 
ing here are observable only in the most slowly rotating stars, 
they reveal that local line profiles are not predicted correctly by 
standard theory, and aside from their interest in allowing us to 
probe the atmospheric velocity fields directly, these effects need 
to be understood (for example through numerical simulations) 
so that they can be correctly incorporated into spectrum syn- 
thesis codes for determining, for example, abundances in stars 
of chemical elements for which only saturated spectral lines are 
available. 

It is observed (as is already well-known) that supergiants in 
this Te range also show clear symptoms of local velocity fields. 
These include large non-zero values of ^, best fit Ve sin ; values 
that increase systematically with equivalent width, and profiles 
of strong lines that are poorly fit by the canonical line model. 
However, such stars do not appear to show the line asymmetry 
found in main sequence A stars, and in the supergiants the line 
profile distortions may be due to g-mode pulsations (Aerts et al. 
'2008) rather than to local convective velocities. 

We conclude by summarising some of the principal observa- 
tional features that theoretical and numerical models should (and 
to some extent do) account for. 
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One major regularity, and one of the strongest pieces of ev- 
idence linking the observed line profile shapes to atmospheric 
convection, is the variation of the strength of the effects dis- 
cussed here with effective temperature. As rises past about 
7000 K, the line profile forms observed change from the solar- 
type shapes to the shapes observed here, with strongly depressed 
blue line wings. As continues to rise, the effect becomes 
steadily weaker until at about ~ 10000 K the line profiles 
become indistinguishable from those predicted by the standard 
model of line shape. As discussed by Landstreet (|1998) . this be- 
haviour is completely consistent with the predictions of mixing- 
length theory. When mixing-length theory is applied to model 
atmospheres appropriate to main sequence A and late B stars, 
it is found that the strength of the super-adiabatic temperature 
gradient in the atmosphere declines rapidly with increasing T^., 
and its driving effect becomes negligible at about the correct 
value of Tf.. Because this behaviour results from a fundamen- 
tal feature of A and B star atmospheres, namely the decline of 
continuum opacity as H becomes more fully ionised with rising 
Je, it is to be expected that it will also be found in numerical 
models. However, because of the high computational expense 
of modelling even a single value, this aspect of the velocity 
field behaviour has not yet been studied numerically; most of the 
modelling has been done for atmospheres with near 8000 K, 
where the effects of the super-adiabatic temperature gradient in 
the lower line-forming region are largest. 

The form of profiles of strong lines, with excess absorp- 
tion far in the blue line wing relative to that in the red wing, 
represents the most severe challenge to numerical modelling at 
present. The observations strongly suggest that the integrated 
form of strong absorption lines is produced by areas of gas rising 
with relatively large velocity, while the descending gas is mov- 
ing at a smaller velocity, which presumably (to conserve mass) 
must cover a larger fractional area. This is the opposite of the so- 
lar case, and is still not found in numerical simulations in spite 
of serious efforts with 2D and 3D codes (e.g. Steffen et al. 120061 
Kochukhovet al. |2007l l. 

A further possible issue arises from the magnitude of the 
"excess" velocities detected, above line broadening due to ther- 
mal motions, pressure broadening, "standard" microturbulence, 
macroscopic broadening by rotation, and finite spectral resolu- 
tion. These may be as large as 10 - 12 kms"' as measured by 
comparing the horizontal discrepancies between observed and 
computed line profiles. Since the speed of sound is roughly 7.5 
to 8.5 km s"' through the line forming region of a main sequence 
A star, the inferred local velocity field may be significantly su- 
personic in places. This does not seem to be a real problem, how- 
ever; Freytag ( 1 995 ) and Kupka et al. ( 2009 ) both find that shocks 
occur in their simulations, indicating clearly the local presence 
of supersonic velocities. 

It is very satisfying that line shape observations are able to 
provide real tests and constraints that challenge efforts to model 
the detailed behaviour of stellar convection. This confrontation 
will certainly help to guide theory eventually to a more physi- 
cally correct description of the phenomenon in stars. 
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